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Rainfall anomalies over southern China are found to be asymmetricly influenced by the Indian Ocean 
Dipole (lOD), with a far stronger influence from positive lOD (pIOD) events. A greater convection anomaly 
and an equivalent-barotropic Rossby wave train response occurs during pIOD events than during negative 
lOD (nIOD) events. Over the Bay of Bengal (BOB) and South China Sea (SCS), an associated low-level 
anomalous anticyclone strengthens the southwesterlies during boreal fall (September, October and 
November, SON), when a pIOD matures. The increased moisture flux gives rise to the anomalously high 
rainfall over southern China. During its developing phase (boreal summer, June, July, and August, J J A), the 
influence of a pIOD event on the contemporaneous rainfall over southern China is weak, but a JJA pIOD 
index is highly correlated with fall rainfall. Therefore, this index can serve as a potential predictor for 
variations of boreal fall rainfall over southern China. 



The Indian Ocean Dipole (lOD) is an air-sea coupled climate mode in the Indian Ocean, characterized by 
opposing sea surface temperature anomalies (SSTAs) in the western and eastern tropical Indian Ocean^"^. It 
not only influences rainfall in the countries around the Indian Ocean, such as Australia, East Africa and 
Indonesia'*"^, but also affects climate variability over remote areas, such as East Asia, Europe and South 
America^"^^. The associated atmospheric teleconnection carries the impact to regions remote from the source 
region^^'^^. The Indian Ocean influences monsoon variability over East Asia via an eastern Indian Ocean/west 
Pacific route, through Rossby wave trains, similar to the Pacific-North American (PNA) patterns^^"^^. 

Many studies have examined the influence of the lOD on the East Asian summer monsoon and China's 
summer rainfalP^. It is suggested that different phases of the lOD mode may change the intensity of the South 
Asian High (SAH) and the western North Pacific anticyclone in summer, and hence variations of monsoon- 
induced rainfall over China^'^^. In contrast, relatively less attention has been paid to the role of the lOD on the 
variations of rainfall over China during boreal fall, when the lOD peaks. A recent study^° showed peak SSTAs for 
the lOD appear in fall and possess a stronger link to the Asian winter monsoon than to the Asian summer 
monsoon. It is suggested that this far-reaching influence during boreal winter is conducted through a wave train 
that is emanated from the Indo-Western Pacific and induced by heating anomalies in this region^\ 

Recent studies^^'^^ have shown that there is an asymmetry between positive and negative lOD events (i.e. pIOD 
and nIOD), with pIOD events reducing rainfall over southern Australia significantly more than a rainfall increase 
from nIOD events. This "impact asymmetry" is consistent with a positive skewness of the lOD^'^^'^^. However, it 
remains unclear how atmospheric circulation over East Asia and rainfall over China respond to the lOD 
amplitude asymmetry. This study examines these issues, and the associated mechanisms. 

Results 

Composite Analysis. The atmospheric circulation over China during warm seasons is characterized by the 
monsoon system^^ Southwesterlies from the Indian Ocean and South China Sea (SCS) advect a large amount 
of moisture to much of China, where rainfall seasonality shows a peak in boreal summer (Fig. la). In boreal fall 
(Fig. Id), an anticyclone flow pattern operates over southern China. The easterly wind to its south and the 
southwesterly wind to its west, converge over southern China, carrying water vapor from the tropical oceans. 
These two moisture transport branches supply the fall rainfall over southern China: one is from the Bay of Bengal 
(BOB) and the other is from the SCS; they originate from the Indian Ocean and the west Pacific, respectively. 
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Figure 1 | China boreal summer (JJA) and fall (SON) rainfall and surface wind vectors, (a), Climatology JJA rainfall (shading) and 850 hPa wind 
(arrows, maximum vector denotes a value of 15 m s~\ and vectors smaller than 1.5 m s~^ are not plotted), (b), Composite JJA rainfall anomalies 
(shading) and 850 hPa wind anomaly (arrows, maximum vector denotes a value of 5 m s~^ and vectors smaller than 0.3 m s~^ are not plotted) for the 
four pure pIOD years (1961, 1967, 1983, and 1994). (c). As in (b) but for the seven pure nIOD years (1958, 1968, 1974, 1980, 1985, 1989, and 1992). (d), 
(e), and (f). As in (a), (b), and (c), respectively, but for SON. White stipples indicate anomalous rainfall is greater than 95% confidence level. All plots were 
produced using licensed Surfer. 



In order to understand the lOD's influence on China's climate, 
composite rainfall and 850-hPa wind anomahes for the pIOD and 
nIOD years are constructed. Based on the classification of a previous 
study^^ the cases of 1961, 1967, 1983 and 1994 in the positive phase 
and those of 1958, 1968, 1974, 1980, 1985, 1989 and 1992 in the 
negative phase are chosen for pure lOD years (without El Nino- 
Southern Oscillation, i.e. ENSO events). The rainfall anomaly com- 
posites of pure pIOD events show that rainfall is above normal over 
southern China (i.e., about south of 28°N), while less rainfall occurs 
over the other parts of China during boreal summer (Fig. lb). In 
contrast, excessive rainfall occurs in most parts of central and south- 
ern China (i.e., about south of 38°N and east of 100°E), while rainfall 
is below normal in northern China (i.e., about north of 38°N) during 
boreal fall (Fig. le). In the pure nIOD years, central China (especially 
the middle reaches of the Yangtze River Valley and the Great Band of 
the Yellow River) experiences a wet summer, but rainfall in southern 
and most parts of northern China declines (Fig. Ic). The fall rainfall 
pattern is almost reversed to a dipole pattern with negative rainfall 
anomalies over southern China but positive anomalies over central 
and northern China (Fig. If). Such a summer lOD-rainfall connec- 
tion is also reported in previous studies^^'^^. 

To further quantify the statistical relation between the lOD and 
China's rainfall, correlation analysis is carried out using the dipole 
mode index (DMI)^ and rainfall anomalies (Fig. 2). Regarding ENSO 
events that are related to the lOD^^'^^, in order to isolate ENSO effects, 
a linear relation of Nino-3.4 with the DMI has been removed from 
the DMI, to obtain a residual index (i.e., DMI | Nino-3.4)- As the lOD is 



weaker in boreal summer compared with boreal fall, no obvious 
correlation between DMI | Nino-3.4 and rainfall anomalies during sum- 
mer is observed (Fig. 2a). However, it still shows that the positive 
phase of the lOD is conducive to a summer rainfall increase in 
southern China while the corresponding negative phase leads to a 
decrease in summer rainfall over the same region, which is consistent 
to the composite analyses (Figs, lb and Ic). As the lOD peaks in 
boreal fall, it has a strong influence on the contemporaneous rainfall 
variations over southern China (Fig. 2c), with correlation coefficients 
significant at the 95% confidence level. 

Although the influence of the lOD in its developing phase on 
summer rainfall over southern China is weak, its influence on fall 
rainfall (i.e., rainfall lags the DMI|Niiio-3.4 by one season) strengthens 
substantially (Fig. 2b). A strong positive relationship between JJA 
DMI|Nifio-3.4 and SON DMI|Nirio-3.4 (Fig- 2d) suggests that the anom- 
alous response of fall rainfall over southern China mainly comes 
from the peaking phase (i.e., SON) of the lOD, because strong JJA 
DMI|Niiio-3.4 tends to be associated with strong SON DMI|Nifio-3.4- 
Consistently, a previous study^° has shown that the matured lOD 
phase has a stronger link to the Asian winter monsoon than to the 
Asian summer monsoon. These results suggest an important point 
that the JJA DMI|Nifio-3.4 can serve as a potential predictor for fall 
rainfall over southern China. 

Asymmetric relationship with SON rainfall and its teleconnec- 
tion. To highlight the asymmetry in the relationship between SON 
rainfall and the opposite phases of the lOD, a scatterplot of SON 
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Figure 2 | Relationship of rainfall with the lOD. (a) , Correlation pattern between JJA rainfall and JJA lOD index (DMI), upon removal of a linear relation 
with the Nino-3.4 index (i.e., DMIlNino-3.4). (b), (c), As in (a) but for SON rainfall and JJA DMI I Nmo-3.4> and for SON rainfall and SON DMIlNino-3.4> 
respectively, (d), Scatterplot of linearly detrended SON DMI I Nino-3.4 against JJA DMI I Nino-3.4- Two sets of linear regression statistics (correlation 
and slope) are shown in (d), using samples with a positive JJA DMI I Nino-3.4 (red), and samples with a negative JJA DMI I Nino-3.4 (blue). Areas within white 
curves in (a), (b), and (c) indicate correlations that are significant at the 95% confidence level, (a), (b), and (c) are plotted using licensed Surfer, (d) is 
plotted using licensed Grapher. 



DMI|Nino.3 4 values versus contemporaneous rainfall anomalies 
averaged over southern China is shown in Fig. 3. For positive 
DMI I Nino-3.4 values, rainfall anomahes tend to be positive 
(although some weak negative values are present), suggesting an 
increase in rainfall, while during nIODs, most rainfall anomalies 
are negative. The influence of a pIOD on rainfall increases is far 
greater than that of a nIOD on rainfall decreases, as indicated by a 
far greater slope during pIOD events. Thus, there is an asymmetry 
between positive and negative events, with a pIOD increasing rainfall 
over southern China significantly more than the rainfall decrease 
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Figure 3 | Response of the China SON rainfall to the SON DMI | Nino-3.4- 
Scatterplot of linearly detrended SON rainfall anomalies averaged over 
southern China (103°-109°E, 24°-28°N) against the SON DMI I Nino-3.4. 
Two sets of linear regression statistics (correlation, slope, and p value) are 
shown, using samples with a positive DMI I Niiio-3.4 (red), and samples with 
a negative DMI I Nino-3.4 (blue). The plot was produced using licensed 
Grapher. 



from the negative phase events, which is probably induced by an 
amplitude asymmetry between the different phases of the lOD (i.e. 
lOD skewness)^"''"-'^ 

What is the relationship between the lOD skewness and the asym- 
metric response of southern China rainfall to the lOD? To examine 
this, we regress grid-point convection and circulation anomalies onto 
the DMI I Nino-3.4- First, we construct their linear relationship using all 
samples (Figs. 4a-4c). Results show a tropically trapped baroclinic 
response to diabatic heating (convective) anomalies over the east/ 
west lOD pole (lODE/IODW), extending northward to the BOB and 
the SCS from the lODE region, to the Arabian Sea from the lODW 
region, respectively (Figs 4a, b). Additionally, in response to the 
convective anomalies over the lOD poles, an eastern Indian Ocean 
and a western Indian Ocean Rossby wave train (Fig. 4a), share a low 
pressure anomaly center stretching from the Tibetan Plateau to the 
northwestern Pacific (Fig. 4b). These Rossby wave-train patterns 
have also been suggested in previous studies^^'^^. The low pressure 
center suggests a weaker SAH than normal. The weaker SAH sup- 
ports a stronger-than-normal summer monsoon over the BOB and 
SCS^^. Consistently, anomalous SON southwesterlies in association 
with the weaker SAH, reach southern China via the BOB and SCS 
(Figs. 5b, 5d), leading to anomalously high rainfall over southern 
China (Fig. 4c). 

Then, we focus on asymmetry with respect to pIOD and nIOD 
events. We use samples with a positive DMI|Nino-3.4j oi" negative 
DMI I Nino-3.4 only in our regression. Circulation anomalies associated 
with pIODs (Figs. 4d and 4e) are better defined than those associated 
with nIODs (Figs. 4g and 4h). This is evident by the presence of 
stronger outgoing longwave radiation (OLR) anomalies over the 
tropics, with corresponding well-defined extratropical wave train 
patterns associated with pIOD events. In fact, much of the wave train 
pattern obtained using all samples (Figs. 4a and 4b) is contributed by 
pIOD events (Figs. 4d and 4e). There is virtually no such signal in 
anomaly patterns associated with nIOD events (Figs. 4g and 4h). This 
asymmetry translates into a direct influence on southern China's fall 
rainfall: the size of positive rainfall anomalies over southern China 
systematically increases with pIOD amplitude (Figs. 3 and 4f). In 
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Figure 4 | Response patterns of circulation to the lOD events, (a), (b), and (c), Circulation anomaly patterns obtained by regressing detrended 
circulation anomalies onto a detrended DMI I Nino-3.4 using all samples averaged over SON, which implies a linear response, (d), (e), and (f), As in (a), (b), 
and (c), respectively, but for samples with positive DMI I Nino-3.4 values (pDMI I Nino-3.4)- (g)> (h), (i), As in (a), (b), and (c), respectively, but for samples 
with negative DMIlNiiio-3.4 values (nDMI I Nino-3.4) - The shading represents the regression coefficients and areas within the thick black contours indicate 
significant correlations at the 95% confidence level, (a), (d), and (g), showthe coefficients associated with outgoing longwave radiation (OLR). In (b), (e), 
and (h), are 200-hPa geopotential height (GH), and in (c), (f), and (i), are rainfall over China. "H" and "L" are marked in (b) to illustrate an eastern 
Indian Ocean and a western Indian Ocean Rossby wave train anomaly centres associated with the lOD. All plots were produced using licensed Surfer. 



contrast, the size of negative rainfall anomalies across the same 
region does not respond to changes in nIOD amplitude in any stat- 
istically significant manner (Figs. 3 and 4i), as the relationship 
between them is not significant at the 95% confident level. 

To further understand the connection of southern China's rainfall 
with lOD events, and to evaluate the associated atmospheric mois- 
ture transportation processes, wind anomaly patterns are con- 
structed by composite analysis (Fig. 5). An asymmetry in wind 
anomalies is also shown, with anomalies far stronger for the pIOD 
samples than for the nIOD years (comparing Figs. 5c-5d to Figs. 5e- 
5f). Consistent with the response to tropical diabatic heating 
described in Figs. 4a and 4d, a strong cyclonic anomaly can be found 
in the upper troposphere over the Tibetan Plateau and southwestern 
China (Figs. 5a and 5c), suggesting a weakening of the upper-level 
tropical easterly winds and anticyclone (i.e. the SAH). On the other 
hand, at 850 hPa, we find an anomalous anti-cyclonic circulation 
extending northeastward from the BOB to the SCS and southern 
China (Figs. 5b and 5d). This is consistent with finding of a previous 
study^°, although their result is based on the data during July to 
September. This circulation anomaly enhances moist southwestly 
winds that blow northward from the BOB and the SCS to southern 
China, facilitating the surplus rainfall over this region. 

Discussion 

We have examined the influence of the lOD on China's rainfall 
during the boreal summer and fall seasons. The results show an 
asymmetry in the impact of pIOD and nIOD events in terms of an 
extratropical response. We find that the extratropical circulation 
response, which is delivered via Rossby wave trains, is greater during 



pIOD events than during nIOD events. During pIOD events, the 
associated low-level anomalous anticyclone found over the BOB 
and the SCS strengthens the southwesterlies. This leads to an 
increased moisture flux that gives rise to the anomalously high pre- 
cipitation over southern China. During nIOD events, the Indian 
Ocean wave train response is far weaker and less well defined. 

The concurrent impact of the JJA DMI | Nino-3.4 on China's summer 
rainfall is weak. However, it has a strong influence on rainfall in the 
ensuing SON season. In fact, the JJA DMI | Nino-3.4 has a comparable 
relationship with SON rainfall to that of the SON DMI | Nino-3.4 with 
SON rainfall. JJA DMI|Niiio-3.4 ^nd SON rainfall possesses a non- 
linear relationship, with a JJA positive DMI|Nifio-3.4 indicating a 
stronger SON rainfall signal than that reflected by a JJA negative 
DMI|Niiio-3.4- In this sense, the JJA DMI | Nino-3.4 serves as a potential 
forecaster for the variations of fall rainfall over southern China. 

Methods 

Data. NCEP/NCAR Reanalysis^° monthly data are used to analyze atmospheric 
circulations that bridge the lOD impacts with climate variability in China. The lOD is 
described using the dipole mode index (DMI)\ defined as the difference of the area- 
averaged SST anomalies in the west lOD pole (lODW, 10°S-10°N, 50°-70°E) and the 
east lOD pole (lODE, 10°S-Eq, 90°- 1 10°E). A monthly in-situ precipitation dataset of 
160 stations in China (Fig. la, white points) from the National Climate Center of the 
China Meteorological Administration is also used to explore the lOD-rainfall 
relationship. All data are in the period of 1951-2012, and monthly anomalies are 
departures from their respective mean values averaged over this period. The monthly 
anomaly data are then linearly detrended and averaged over the seasonal months to 
obtain the seasonal anomaly data. 

Analysis technique. The basic analysis technique used in this study is linear 
regression (or equivalently, correlation). Before regression analysis, in order to 
exclude ENSO effects, a linear relation of Nino-3.4 with the DMI is removed from the 
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Figure 5 | Boreal fall (SON) 200-hPa and 850-hPa winds during pure lOD years, (a), Composite SON 200-hPa wind anomalies (arrows, maximum 
vector denotes a value of 6.0 m s~^ and vectors larger than or equal to 0.5 m s~^ are plotted in blue while anomalies smaller than 0.5 m s~4n black) for the 
11 pure lOD years (i.e., including the four pure pIOD years (1961, 1967, 1983, and 1994), and the seven pure nIOD years (1958, 1968, 1974, 1980, 
1985, 1989, and 1992)). (c), (e). As in (a), but for the pure pIOD years and the pure nIOD years, respectively, (b), (d), and (f), Asin (a), (c), and (e), but for 
850-hPa winds, respectively. All plots were produced using licensed Surfer. 



DMI to obtain a residual index (i.e., DMI|Nino-3.4)- Regression coefficients (and the 
corresponding correlations) are first calculated by regressing anomalous circulation 
fields onto the DMI|Nino.3 4 using all samples, assuming a symmetric response with 
respect to positive and negative phases. Regression coefficients are then calculated 
using samples containing only pIOD or only nIOD values. Statistical significance of 
the regression anomalies is estimated based on significance of the associated 
correlation coefficients, which is when p value < 0.05. 

Seasonality. All calculations in this paper are based on data in boreal summer (June, 
July, and August, or J J A) and in boreal fall (September, October, and November, or 
SON), when lOD events develop and peak respectively. 

Graphic software. All maps and plots were produced using licensed Surfer or 
Grapher. 
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